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ON 

A comparison is made between two variable flavor number schemes which 
describe charm quark production in deep inelastic electron-proton scattering. 
JL . In these schemes the coefficient functions are derived from mass factorization 

of the heavy quark coefficient functions presented in a fixed flavor number 
scheme. Since the coefficient functions in the variable flavor number schemes 
have to be finite in the limit m^Owe have defined a prescription for those 

^ processes where the virtual photon is attached to a light quark. Furthermore 

one has to construct a parton density set with four active flavors (u,d,s,c) out 
of a set which only contains three light flavors (u,d,s). In order a 2 s the two 
sets are discontinuous at fi = m c which follows from mass factorization of 
the heavy quark coefficient functions. The charm component of the structure 
function i*2,c is insensitive to the different variable flavor number schemes. 
In particular in the threshold region they both agree with the description in 
fixed order perturbation theory presented in a three flavor scheme. However 
one version does not lead to a correct description of the threshold behavior 
of the longitudinal structure function F L)C . This happens when one requires 
a non-vanishing zeroth order longitudinal coefficient function. 

PACS numbers: 11.10 Jj, 12.38Bx, 13.6011b, 13.87Ce. 



1 Introduction 

Charm quark production is one of the important reactions used to extract 
the gluon density f g (x,fi 2 ) of the proton in deep inelastic lepton-hadron 
scattering, especially when the Bjorken scaling variable x is small. However 
this is only true when the deep inelastic process is of the neutral current 
type and the charm component of the proton wave function is negligible. In 
this case the charm quark is produced in the so-called extrinsic way. For 
neutral current processes with only light partons in the initial state this 
means that the Born approximation in perturbative QCD is given by the 
virtual vector-boson gluon-fusion process p]. Notice that the light partons 
consist of the gluon and the three light flavors u, d, s together with their 
anti-particles. Furthermore if the virtuality of the exchanged vector boson 
in deep inelastic lepton-hadron scattering satisfies Q 2 <C M§ then the vector 
boson is represented by the photon only and the contribution of the Z-boson 
is negligible. Extrinsic charm production also receives next-to-leading order 
(NLO) contributions from boson-quark subprocesses, which could hamper 
the extraction of the gluon density. Fortunately this is not the case at HERA, 
where the experiments 0, are carried out at small x, because the gluon 
density overwhelms the light flavor densities completely. Moreover the NLO 
quark initiated processes are suppressed by at least one power of the strong 
coupling constant a s (/x 2 ) with respect to the Born contribution to the boson- 
gluon fusion reaction. The quantity /i in the running coupling constant and 
the parton densities represents both the renormalization and factorization 
scales respectively, because it is convenient to chose them to be equal. 

In the literature one has adopted two different treatments of extrinsic 
charm production, which are known as the massive and massless charm de- 
scriptions. The former, advocated in Q, treats the charm quark as a heavy 
quark (with mass m c ) and the cross sections or coefficient functions have 
to be described by fixed order perturbation theory. Notice that due to the 
work in [Rj the perturbation series is now known up to second order and 
the NLO massive charm approach agrees with the recent data in || and [[|. 
The latter treatment, which has been rather popular among groups which 
fitted parton densities to experimental data, treats the charm quark as a 
massless quark so that it can be represented by a parton density / c (x,/i 2 ), 
with the boundary condition f c (x,fi 2 ) = for // < m c . Although at first 
sight these approaches are completely different they are actually intimately 



related. It was shown in || that the large logarithms of the type ln(Q 2 /m^), 
which appear in the perturbation series when Q 2 3> m 2 c) can be resummed 
in all orders. The upshot of this procedure is that the charm components of 
the deep inelastic structure functions F ijC (x,Q 2 ,m 2 ), where i = 2,L, which 
in the first approach are written as convolutions of heavy quark coefficient 
functions with light parton densities, become, after resummation, convolu- 
tions of light parton coefficient functions with light parton densities which 
also include a charm quark density. This procedure leads to the so-called 
zero mass variable flavor number scheme (ZM-VFNS) for Fi tC (x, Q 2 ) where 
the mass of the charm quark is absorbed into the new four flavor densities. 
To implement this scheme one has to be careful to use quantities which are 
collinearly finite in the limit m c — > 0. From the above considerations it is 
clear that the first approach is better when the charm quark pair is produced 
near threshold because the mass of the quark is important in this region and 
it cannot be neglected. On the other hand far away from threshold, where 
also Q 2 ^> m 2 , the large logarithms above dominate the structure functions 
so that the second approach should be more appropriate. Both approaches 
are characterized by the number of active flavors involved in the description 
of the parton densities which are given by three and four respectively. There- 
fore one can also speak of three and four flavor number schemes (TFNS and 
FFNS respectively). Each scheme has a different gluon density so that the 
momentum sum rule is always satisfied. 

As most of the experimental data occur in the kinematical regime which 
is between the threshold and the region of large Q 2 a third approach has been 
introduced to describe the charm components of the structure functions. This 
is called the variable flavor number scheme (VFNS). A first discussion was 
given by Aivasis, Collins, Olness and Tung |7|], where a VFNS prescription 
called ACOT was given in lowest order only. The ACOT results were com- 
pared with the NLO results in || . We will give our NLO version of a VFNS 
scheme in this paper and we call it the CSN scheme to distinguish it. A dif- 
ferent approach, generalized to all orders, was given in the papers by Buza, 
Matiounine, Smith and van Neerven [|J,||, which we denote by BMSN. Fi- 
nally another version of a VFNS for the charm component of the structure 
function was presented by Thorne and Roberts in [Hj, which will be called 



the TR scheme. Note that a proof of factorization to all orders for the total 
structure function, which includes charm and light parton production, was 



recently given in [11 



The difference between the various versions can be attributed to two in- 
gredients entering the construction of a VFNS. The first one is the mass 
factorization procedure carried out before the large logarithms can be re- 
summed. The second one is the matching condition imposed on the charm 
quark density, which has to vanish in the threshold region of the production 
process. It will be one of our goals to elucidate these differences in the next 
Section. Another problem, which was not clarified in the papers above is 
that the mass factorization cannot be carried out on the level of the charm 
components of the structure functions alone, because one also needs contri- 
butions coming from the light parton components of the structure functions. 
The latter can be attributed to all heavy charm quark loop contributions 
to gluon self energies, which appear in the virtual corrections to the light 
parton coefficient functions. These corrections have to be combined with 
contributions from gluon splitting into heavy charm anti-charm quark pairs, 
which belong to the charm components (not the light quark components) 
of the structure functions. In this paper we will give a much more careful 
analysis than has been done previously in the literature. Another aspect of 
any VFNS approach is that one needs two sets of parton densities. One set 
only contains densities in a three flavor number scheme whereas the second 
one, which also includes a charm quark density, is parametrized in a four 
flavor number scheme. Both parameterizations have to satisfy the relations 
quoted in ||. At this moment the latter set is not available in the literature 
and we would like to fill in this gap. Starting from a three flavor number set 



of parton densities recently published in [12| we will construct a four flavor 
number set of densities satisfying the relations in || . 

In Sec. II we give a general discussion of the CSN description for heavy 
quark electroproduction, and explain the problems with mass factorization, 
collinear singularities and threshold dependence in the heavy flavor compo- 
nents of the structure functions. We then specialize to charm quark elec- 
troproduction in Sec. Ill, working to second order in the running coupling 
constant a s (/z 2 ). We first present details about the charm quark density. 
Next numerical results are shown for the structure functions in the various 
schemes. Analytic results for the contributions from the Compton scattering 
reaction with an invariant mass cut are relegated to an Appendix. Finally 
we want to emphasize that we only consider inclusive charm quark produc- 
tion in this paper. Exclusive charm production which involves transverse 
momentum and rapidity distributions will be dealt with in another paper. 



2 Discussion of variable flavor number schemes 

In this section we discuss two different representations of the deep inelastic 
structure functions in variable flavor number schemes. One is proposed here 
(CSN). The other (BMSN) was proposed in || and ||. The former starts 
from mass factorization of the exact heavy quark coefficient functions whereas 
the latter only applies this procedure to the asymptotic expressions for these 
functions. In both schemes the special role of the heavy quark loop contribu- 
tions to the light quark coefficient functions in combination with heavy quark 
production via gluon splitting was overlooked. This will be repaired in this 
paper. Furthermore in both schemes there is a lot of freedom in the choice 
of matching conditions, which are needed to connect the structure functions 
presented for rif and rif + 1 light flavors. Different matching conditions lead 
to different threshold behaviors, which have consequences for the description 
of the structure functions at small Q 2 and large x. 

Limiting ourselves to electroproduction, where deep inelastic lepton-hadron 
scattering is only mediated by a photon, the light parton components of the 
structure functions are defined by 
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where ® denotes the convolution symbol in the parton Bjorken scaling vari- 
able z. In this expression the C^ (i — 2, L; k — q, g) denote the light par- 
ton coefficient functions and the e& represent the charges of the light flavor 
quarks. The quantities C^ RT only contain the heavy quark loop contribu- 
tions to the light parton coefficient functions. Furthermore fg(nf, fj 2 ) stands 
for the gluon density while the singlet (S) and non-singlet (NS) light quark 



densities, with respect to the SU(n/) flavor group, are defined by 
fk+k( n f, V?) = fk(n f , //) + fk(n f , /x 2 ) 

n f 

fq( n f^ 2 ) = Yl fk+k(n f ,fl 2 ) 

k=l 

/f(n/,/i 2 ) = /*+*("/, /i 2 ) ~ ^f!(n f ^ 2 ) ■ (2.2) 

Finally we have set the factorization scale equal to the renormalization scale 
/i. The light parton coefficient functions have been calculated up to order a 
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in [ 13(1 . The contributions to C 4 V fc RT appear for the first time in second order 
perturbation theory and can be found in [13]. For our further discussion it 
will be convenient to distinguish between the numbers of external and internal 
flavors. The former refers to the number of light flavor densities whereas the 
latter denotes the number of light flavors in the quark loop contributions 
to the virtual corrections. They are not necessarily equal. Some of the 
coefficient functions have the external flavor number as an overall factor. To 
explicitly cancel this factor we have defined the quark and gluon coefficient 
functions in Eq. ([2.1|) as follows 
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C S (-/> $) = n f C?*(n f , ^) , Cl(n f , ^) = n f Cf !g (n f , ^) , (2.3) 
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where PS represents the purely singlet component. Hence the remaining rif 
in the argument of the coefficient functions marked with a tilde denotes the 
number of internal flavors. The same holds for the n/ in the parton densities. 
However the argument nf in the structure functions is external and it refers 
to the number of parton densities appearing in their expressions. The parton 
densities satisfy the renormalization group equations. If we define 

D = fi— + (3{n f ,g)—, g = g{n f ,fi 2 ) (2.4) 



then 
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Figure 1: Lowest-order photon-gluon fusion process 7* + g —* Q + Q con 
tributing to the coefficient functions H^ . 
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where 7^ represent the anomalous dimensions of the operators in the operator 
product expansion (OPE). 

The heavy flavor components (Q = c, b, t,Q = c, b,t) of the structure 
functions F% and Fl arise from Feynman graphs with heavy flavors (Q and 
Q with mass m) in the final state and are given by 
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where eg represents the charge of the heavy quark. Furthermore L it k and H^ 
(i = 2, L; k = q, g) represent the heavy-quark coefficient functions which are 
exactly calculated order by order in perturbation theory. In Figs. |TJ-^| we have 
shown some of the Feynman diagrams contributing to the coefficient functions 
up to order a 2 . Like in the case of the light-parton coefficient functions C^k 
they can be split into (purely)-singlet and non-singlet parts. The distinction 
between L^ and H^k can be traced back to the different (virtual) photon- 
parton heavy-quark production mechanisms from which they originate. The 




Figure 2: Virtual gluon corrections to the process 7*+g — > Q+Q contributing 
to the coefficient functions H, '} . 



functions L^k, Hi,k are attributed to the reactions where the virtual photon 
couples to the light quarks and the heavy quark respectively. Hence L^ 
and H ik in Eq. fl2.6| ) are multiplied by e 2 . and eg respectively As has 
been mentioned in the introduction the heavy quark coefficient functions 
contain large logarithms of the type ln l (Q 2 /m 2 ) when Q 2 3> m? which can 
be removed from the former by using mass factorization. To do this we first 
have to split the heavy quark coefficient functions L^ into soft and hard 
parts 
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where A is a cut on the invariant mass Sqq of the heavy quark pair. The 
cut can be determined by experiment. It is chosen such that in the limit 
all mass singularities reside in the soft parts so that the hard parts 



m 



are collinearly finite. Taking Fig. |5| as an example we mean by hard that 
one detects a QQ-pair with a large invariant mass which is experimentally 
observable if Sqq > A. In the case Sqq < A the QQ-pair is soft and becomes 
indistinguishable from other light parton final states which contain contribu- 
tions from virtual heavy quark loops. Next we add the soft parts to the other 
contributions to i^ LIGHT in Eq. (|2.1|) and the mass factorization proceeds like 
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Here Cj,;^ are those parts of the light parton coefficient functions Cij which 
contain the heavy quark loops. The hard parts of L^k are left in F®q ACT i n 
Eq. (|2.6| ) and do not need any mass factorization. Furthermore we have the 
condition that the dependence on the parameter A cancels in the sums so 

f>CSN/_ Q Q \ _ ^CSN.SOFT/ A Q Q \ 

+ Lgf»(n / ,A,^ > 5), (2.9) 

where /i in the hard parts only represents the renormalization scale. The 
coefficient functions H^ satisfy the relations [] 

H hk {n f ^%) = A lk {n f ^)®C^{n f ^%) k, I = Q, q,g .(2.10) 

Notice that mass factorization applied to the functions H iq and H ig occurring 
in Ffq ACT ( |2.6|) leads to the coefficient functions C^q N - The latter also follow 
from mass factorization of the functions H it q which represent processes with 
a heavy quark in the initial state. The quantities H it Q, which do not appear in 
Ffn :ACT ; together with the corresponding operator matrix elements (OME's) 
Aqq are characteristic of variable flavor number schemes. The procedure 
above transfers the logarithms ln 4 (/z 2 /m 2 ), appearing in Lf® FT and Hi^, to 
the heavy quark operator matrix elements A^q and AQ k . The latter are 
defined by (see |§ for details of renormalization and mass factorization) 

Ak,Q = (k\Oi(0)\k) k,l = q,g, 

A Qk = (k\O Q (0)\k) k = Q } q,g. (2.11) 

Note that the 0\ are the light quark and gluon operators and in Ad-^q we only 
retain contributions from subgraphs which contain heavy quark (Q) loops. 

1 In order to help the reader, who is more familiar with the notation in H], one can 
make the following comparison. For instance Eq. (7) in the latter reference fg is equal 
to our Ag>. Similarly in Eq. (8) Ea^? = H S and £a<4q = H l% 
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Figure 3: The bremsstrahlung process l* + g-^Q + Q + g contributing to 
the coefficient functions H^g . 



The quantity Oq represents the heavy quark operator. Here we want to stress 
that the operators are sandwiched between quark and gluon states. This will 
cause mass singularities of the type ln*(/x 2 /m 2 ) to appear in a similar way as 
they appear in partonic cross sections. 

The heavy quark coefficient functions defined in the CSN scheme in Eqs. 
( |2.8|) - (|2.10|) are collinearly finite and tend asymptotically to the massless 
parton coefficient functions presented in Eq. (|2.1|) i.e. 
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In particular we have 



C*,*(n/, %) + lim C^(n f} % JL) = C i>k {n f + 1,%), (2.13) 



so that the number of internal flavors is enhanced by one unit. 

The CSN scheme above has similarities with the VFNS schemes proposed 

in i, 
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The decomposition of the L^k into soft and hard parts has not 
been discussed previously. However one must address this issue because the 
mass singularities in L^*. and CYl RT separately have such high powers in 
\n(Q 2 /m 2 ) that they cannot be removed via mass factorization. Another 
feature is that in the limit m — > the final state invariant energies in the 
reactions which contribute to these two types of coefficient functions become 
equal. Hence L iy k and C 4 V fc RT have to be added so that the leading singular- 
ities cancel and the remaining ones are then removed by mass factorization 
according to Eq. (|2.8|). Notice that the the total coefficient functions L ik 



11 



cannot be moved to i^ LIGHT because this would contradict the definitions of 
the latter structure functions where only light partons are observed in the fi- 
nal states. Therefore it is sufficient to transfer the Lf ° FT to the i^ LIGHT since 
they contain the same mass singularities as the L^. If A is chosen small 
enough, the heavy quarks are unobservable in a measurement of _F i LIGHT . 
Note that the problem of the separation of L^ into soft and hard parts is 
not needed for the total structure function F LIGHT -+- _p E ^ ACT . The all order 



mass factorization of the latter is shown in [|Tl|] . 

To illustrate the procedure above we carry it out up to order a 2 . The 
coefficients in the series expansion are defined as follows 
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Up to second order the mass factorization relations become 

}\ .SOFT,NS,(2)/ A &_ &_\ _ ,NS,(2)/_/i_ 2 
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with 
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7 J TOO *• TOO *■ ■* 



Here F^ 2 \Q 2 /m 2 ) denotes the two-loop vertex correction in Fig. |6|. This 
function satisfies the decoupling theorem which implies that it vanishes in 
the limit m — > 00. The heavy quark coefficient functions L ik have been 



calculated in [15] and, after their convolution with the partonic densities, 
yield contributions to the structure functions which behave asymptotically 
like In (Q 2 /m 2 ). These logarithms are canceled after adding F^(Q 2 /m 2 ) 
in Ref. |14j] to L ik ' which contains the same mass singularities as 

L ik A (see the remark below Eq. (|2.7| )). To obtain the hard and soft parts 
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Figure 4: Bethe-Heitler process 7* + q(q) — > Q + Q + g(g) contributing to 

PS (2) 

the coefficient functions H iq . The light quarks q and the heavy quarks Q 
are indicated by dashed and solid lines respectively. 



we divide the integral over sqq in the graphs of Fig. |5] in two regions i.e. 
s > Sqq > A and A > Sqq > 4 m 2 which we denote by HARD and SOFT 
respectively. Here Sqq and s denote the CM energies squared of the QQ 
system and the incoming photon-parton state respectively. The hard and 
soft parts are presented in the Appendix. Finally /3o,q = —2/3 denotes the 
heavy quark contribution to the lowest order coefficient of the /3-function in 
Eq. (2.4). We must change the running coupling constant when we change 
schemes. 

The mass factorization of the heavy quark coefficient functions H^ is 
simpler. Here we get 

H PS,(2)^Q^_ Q^ = ^PS,(2)^ 2 ^CSN,NS,(0)^Q 2 >| | g CSN,PS,(2) ( Q 2 Q 2 \ 

(2.17) 
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(2.19) 



13 



P 






P 



q q 

Figure 5: Compton process ■~f*(q) + q(p) — > Q(Pi)+Q(P2)+q(p') contributing 
to the coefficient functions L„-„ . The light quarks q and the heavy quarks 



hi 



Q are indicated by dashed and solid lines respectively (s = (p + q) 2 , s 
(pi +P2) 2 see text). 
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In the expressions above we have only given the arguments on which the 
coefficient functions and operator matrix elements depend, like n/, Q 2 /m 2 , 
or Q 2 /fi 2 (at least up to that order in perturbation theory). Furthermore 
we have dropped the convolution symbol when the corresponding coefficient 
function behaves as a 5- function of the type 5(1 — z). The heavy quark 
coefficient functions correspond to the following processes 
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1* + 9^Q + Q Fig. 1 
l* + 9^Q + Q + 9 Figs. 0| 
q(q) ^Q + Q + q(q) Fig. | 
q{q) ^Q + Q + q(q) Fig. | 
7* + Q -+ Q 
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Figure 6: Two- loop vertex correction to the process 7* + q — > q con- 
taining a heavy quark (Q) loop. It contributes to C iq ' (Q 2 /m 2 ) = 

F^{Q 2 /m 2 ) C\ 
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i,q ■ 
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NS,(1) 
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Q^Q + g Fig. 0. 



(2.21) 



In the reactions above the virtual corrections to the lowest order processes are 
implicitly understood. The coefficient functions L i 



^' (2) , flJf'W and HF< 



i,g 



whereas the H, n are 



computed in the MS'-scheme, can be found in 
computed in the context of QED in jRJ. The MS-scheme is also chosen 
for the OME's in Eqs. ( |2.15 )-( 2.20 ) which are computed up to order a 2 in 
U and [15], This also holds for Aqq in Eq. ( |2.18|) which is presented 
in the context of QED in p7]. Furthermore the running coupling constant 
appearing in the quantities above contains nf active flavors. 

The mass singular logarithms of the type ln(/i 2 /m 2 ), appearing in the 
OME's above, are absorbed by the light parton densities. This procedure 
leads to parton densities which are represented in the rif + 1 light flavor 
scheme. For the light parton densities one obtains 
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The parton density representing the heavy quark in the iif + 1 flavor scheme 
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Finally the gluon density in the nf + 1 flavor scheme is 



r 



fg(n f + l,fi 2 ) = A b gqjQ (n f ,^)®f»(n f ,fi 2 ) 



m* 



// 



m 



(2.23) 



+ A 99 d n f^)®f>f^)- ( 2 - 24 ) 



One can check (see || ) that the new parton densities satisfy the renormal- 
ization group equations in Eq. Q2.5|) wherein all quantities nf are replaced 
by nf + 1. Up to order a 2 the above relations become 
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/>/ + l^ 2 ) = /!("/, A* 2 ) +a.(n/,A« 3 )4J3(^) ® />/,** 2 ) 



+a 2 (n /; /i 2 ) [^S(^) ® />/, Z^ 2 ) + A S g i 

lib 



WrJL 

in 



(>i — i&f%(nf,H 2 ) 



(2.27) 
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Figure 7: Order a s corrections to the process 7* + Q — > Q and the reaction 
7* + Q — »■ Q + g contributing to the coefficient functions H { q' ( ' . 



Notice that in passing from an n/-flavor to an n/+l-flavor scheme the running 
coupling constant becomes 

r v 2 1 

a 8 (nf + l,n 2 ) = a s (n f ,n 2 )\l - a s (n /; /i 2 )/3o,Qln( — )J , (2.28) 

fit 

which has to be used in all expressions for the structure functions in the 
CSN scheme. Using the mass factorization relations in Eqs. ( |2.8| ), ( 2.10| ) and 
the redefinitions of the parton densities in Eqs. ( |2.22|) -(2.24) we obtain from 
Eq. (|2.6| ) the heavy quark components of the structure functions in the CSN 
scheme are 



itT(n / + l,A,Q 2 



m 



Q 



Q 2 Q' 



t I 1 1 2\ ^ I^CSN.NS/ , 1 V V 

fq+Q^f + 1, /i ) ® C i>Q (n, + !,—,_ 



mr /r 



, ^CSN,PS/ . 1 Q Q \\ 



1 V t l i 1 2\ „ ^CSN.PS/ . -, Q Q 



1=1 



i -fS/ ,1 2\ ,«, r>CSN,S/ , -t Q Q 



+ E^ 



fc=l U=l 



Ef / . -, 2\ ^ r HARD, PS/ * Q Q 



77r /i^ 
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i fSf ,1 2\ ~ r HARD,S/ A Q Q 

+f g (n f + l,ii)®L i)g (n f ,A,—,— 



. p i . -, 2\ ~ r HARD,NS 



/ A Q 2 Q 2 



(2.29) 



In a similar way one obtains from Eq. ( |2.1|) the light parton components of 
the structure functions in the CSN scheme are 



^CSN.LIGHT 



k=i 



n f 



E /i+K"/ + 1 



i=i 



2\ -o (nP$( Q \ i A>CSN,SOFT,PS / A Q Q \ 1 



-f*{n f + l,if 



+A+fe(^/ + l^ 2 



C l 9 ( ra /'-^ 



, ^CSN.SOFT.S/ a Q Q A 



??? 






+ c. 



CSN,SOFT,NS 
i,q,Q 



(n f ,A, 






(2.30) 



Up to a given order Eqs. (|2.29|) and ( |2.30|) do not differ from the structure 
functions presented in Eqs. (|2.6|) and ( |2.1D respectively as long as one uses 
fixed order perturbation theory. This can be checked up to second order 
when the coefficient functions in Eqs. (|2.29| ) and ( |2.30|) are substituted using 
the mass factorization relations in ( |2.15|) - (|2.20|) . The difference arises if the 
logarithms of the type ln 4 (/z 2 /m 2 ), which show up in the parton densities, are 
resummed using the renormalization group equations in Eq. ( |2.5| ). This re- 
summation induces corrections beyond fixed order perturbation theory which 
become noticeable for fi 2 ^> m 2 . On the other hand we do not want that 
the resummation bedevils the threshold behavior of the structure functions. 
In this region the best representation is still given by Eqs. ( |2.1| ) and (|2.6|) . 
Therefore one has to look for a scale at which expressions ( |2.29| ) and ( 2.30Q 
coincide with those given by fixed order perturbation theory in (|2.6|) and 
( |2.1| ) respectively. Finding this scale is the most important issue in CSN 
as we will show below. Both expressions FpQ N and F t ' are renor- 

malization group invariants. Hence they satisfy the renormalization group 
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equations (see Eq. (pT 



D iff = , D ^SN.LIGHT = Q (2 31) 

The same holds for the total structure function in the variable flavor number 
scheme which is defined as 

if Sr V/ + 1, Q 2 , m 2 ) = if SN ' LIGHT (n /5 A, Q 2 , m 2 ) 

+F^(n f + l,A,Q\m 2 ). (2.32) 

One can now show that for large Q 2 , if (n/+l, Q 2 , m 2 ) turns into the same 
expression as Eq. ( |2.ip where Uf in the light parton coefficient functions Cij. 
is replaced by rif + 1 and C 4 V fc RT = f° r the rif + 1 heavy flavor piece. 

After having discussed the general procedure to construct CSN scheme 
structure functions we now turn to the practical issues. For asymptotic values 
of Q 2 , far above the heavy QQ threshold at (1 — x)Q 2 /x = 4m 2 , all coefficient 
functions Cff N in Eq. ( |2.29| ) can be replaced by the light parton coefficient 
functions so that, after having removed L^ RD , one gets the heavy quark 
components of the structure functions in the so-called zero mass variable 
flavor number scheme (ZM-CSN). However near threshold at low Q 2 and 
large x there is a problem, which has not been solved satisfactorily in the 
literature. In this region one would like the ifg N to vanish in the same way 
that the FpQ ACT vanish. Unfortunately the coefficient functions (f q N do not 
vanish in the threshold region due to the presence of the OME's Aq k and 
the functions Hfg which describe processes with ONE heavy quark in the 
final state (see Eq. ( |2.21| ) ) contrary to Hi^ g and Hi q which originate from 



reactions with at least TWO heavy quarks (Q and Q) in the final state. Only 
the latter functions have the correct threshold behavior. In the literature two 
ways have been proposed to obtain reasonable threshold behaviors. The first 
one was given in a paper by Aivasis, Collins, Olness and Tung [0], which will 
be denoted as the ACOT boundary conditions. The second one was proposed 
in a paper by Thorne and Roberts [|IT| , which we shall call the TR boundary 
conditions. In both approaches the structure functions have the properties 

F hQ (Q 2 ,m 2 )=F% ACT (n f ,Q 2 ,m 2 ) for Q 2 < m 2 , 
F hQ (Q 2 ,m 2 )=F™ N (n f + l,Q 2 } m 2 ) for Q 2 > m 2 , (2.33) 
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where the parton densities satisfy the boundary conditions at /i 2 = m 2 
fk+i(n f + l,x,n 2 ) = f k+ k(n f ,x,fi 2 ), 

f g (n f + l,x,n 2 ) = f g (n f ,x,/j, 2 ). (2.34) 

Notice that there is no relation between the scale /x, chosen in these bound- 
ary conditions, and the production threshold of heavy quarks (1 — x) Q 2 /x > 
4 m 2 . If one e.g. takes /i 2 = Q 2 and Q 2 < m 2 all terms in Eq. ( |2.29| ), where 



the heavy quark density is multiplied with Cj C q N , vanish in spite of the fact 
that heavy flavors are still produced as long as x < Q 2 /(Q 2 + 4 m 2 ). On the 
other hand it is possible that Q 2 > m 2 and x > Q 2 /(Q 2 + 4 m 2 ) which implies 
a non-vanishing heavy quark density with no heavy quark pair production. 
Therefore the value for Q 2 chosen in Eq. ( |2.33j ) is a little arbitrary and no 



condition is imposed on x. Since we do not have any alternative we shall 
choose the same value of Q 2 as in Eq. (|2.33j ) above which F®q ACT turns into 
FpQ N . Another problem is is that the boundary conditions in Eq. (|2.34|) do 



not agree with the relations in Eqs. ( |2.22| )-( f2.24| ) if the computations are 
extended beyond order a 2 which happens in the MS-scheme. In this case 
the continuity at /i 2 = m 2 is changed into a discontinuity. Finally there is 
a problem with the heavy quark density if we choose an arbitrary value for 
/i 2 in /q+q(^/ + 1, x, /i 2 ). In deep inelastic scattering one very often chooses 
/i 2 = Q 2 . If Q 2 < m 2 one has to know /g + g(n/ + 1, x, /i 2 ) for values /i 2 < m 2 
which are not specified in Eq. ( |2.34| ). Furthermore the OME's do not vanish 
for /z 2 < m 2 so that Cf^ ^ H i>k . (see Eqs. (|2~il)-(|2~20|)). Hence in Eq. 
(|2~3l ) F i:Q ^ Fffi ACT for Q 2 < m 2 so that the threshold behavior will be 
spoiled. Therefore one has to avoid chosing /i 2 < m 2 which can be achieved 
by the prescription given by ACOT in 



2 2 , ; r\i ' m 



2\ n 



^ = m ' + kQ z 11 - — \ for Q z >m 2 , 

fi 2 = m 2 for Q 2 <m 2 , (2.35) 



^^Note that in the existing parton density sets fn+ninf + l,x,[x) also vanishes for 



[l 1 < to 2 . 
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with k = 1/2 and n = 2. In this way one gets Cp^ = Hi^ for Q 2 < m 2 
at least up to order a s . For higher orders one has to use the relations in 
Eqs. ( |2.22p -( [2~2~iD instead of those given in Eq. fl2.34| ) as the latter only hold 



up to order a s . The new conditions are presented up to order a 2 in Eqs. 
(pH)-(2.27). 

In the TR prescription one first chooses fi 2 = Q 2 and then requires 

F^(n f + 1, Q 2 , m 2 ) \ Q , =m , = if XACT (n /; Q 2 , m 2 ) \ Q2=m , , 

dF^{n f + l,Q 2 ,m 2 ) dif XACT (n /5 Q 2 ,m 2 ) 



d ln(Q 2 /m 2 ) lv ~ m d ln(Q 2 /m 2 ) 



\Q 2 =m 2 ■ 

(2.36) 



Using the mass factorization relations in Eqs. (|2.17|) - (|2.20|) the TR boundary 
conditions lead to new heavy quark coefficient functions C^q which have 
nothing to do with the reactions in Eq. (|2.21|) . For instance the lowest order 



coefficient functions in the TR prescription corresponding to the reaction 7* + 



Q —* Q in Eq. (|2.21|) vanish at the QQ threshold although this process only 
contains one heavy quark in the final state. Moreover, as already admitted 
by the authors in pTOfl , this procedure breaks down beyond order a s because 



there are more coefficient functions than relations between them. Another 
problem is that it is unclear in which subtraction scheme one is working 



since the subtraction terms in |L0| have nothing to do with the usual OME's 
except in the limit Q 2 3> m 2 . For example in order a s the subtraction term 
needed for H { q in Eq. (|2.18|) is not given by Aqq . The same holds for 

H i g in Eq. ( |2.19| ), which gets another subtraction than given by Aq 



From the theoretical viewpoint the boundary conditions in Eq. ( |2.36|) seem 
very unattractive to us because the relations between the coefficient functions 
^ C fc N (^ = Qi Qi d) anc i the actual parton reactions are broken and the scheme 
is unknown. Notice that the parton densities in [|TI]] are still presented in the 
MS-scheme. 

A different VFNS from that discussed above has been proposed by Buza, 
Matiounine, Smith and van Neerven in ||, 0, which we call the BMSN 
scheme. In the latter it was advocated that only when the large logarithms 
dominate the heavy quark coefficient functions do they have to be removed 
via mass factorization so that they can be resummed via the renormalization 
group equations. In the BMSN scheme we need the asymptotic heavy quark 
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coefficient functions defined by 

12 r& 



lim H hk {n f , — , — ) = H t ° (n f , —, — ) , 

Q 2 ^>m 2 m z fi z m z \i z 



(2.37) 



which behave like 

H ASVMP,,„ ( , Vi ? Q!) „ „, ^ ^ , n „ (?) , n , g) (2 3g) 

'" r n+j<l " i r 

with a similar behavior for Lff YMP . In the BMSN scheme Ffq ACT is given by 
Eq. ( |2.6| ) except that L^ — ► -^ffc" 110 analogous to the CSN. The asymptotic 
heavy quark structure functions, denoted by F^ YMP , are given by the same 
expressions as presented for F^q ACT where now all exact heavy quark func- 
tions are replaced by their asymptotic ones. Up to second order the latter 

u t A ■ [FBI rr, p ,■ r SOFT,ASYMP,NS,(2) , T HARD,ASYMP.NS,(2) 

can be found m |p.q| . 1 he functions L iq and L i 

are given in Appendix A. In the BMSN scheme the charm components of the 

structure functions are defined as 



F$™(n f + 1, A, Q\ m 2 ) = F^(n h A, Q 2 , m 2 ) 
-KT MP (n f , A,Q 2 ,m 2 ) + F™ F (n f + 1, A,Q 2 ,m 2 ) 



(2.39) 



with 



2 ™2\ 



F™*(n f + l,A,Q\m 



Q 2 



/ g V/ + l^ 2 )®^?(n/ + l,^)+jf(n/ + l,A« 2 ) 



/'• 



Q 



Cl (n f + l,^) + f Q+Q (n f + 1, /i 2 ) ® C™ (n, + 1 



Q s 



+ E 



*!=i 



HARD,ASYMP,PS 



.1=1 

, f%( ,i 2\ r HARD,ASYMP,S/ a 

+ fg( n f + 1 ^ )® L i,9 (%> A > 

. t i , -, 2\ Q r HARD,ASYMP,NS/ a 



/ A Q 2 Q 2 



2 /02 



g^ g_ 

9 ' 9 
771 /i z 



2 /n2 



cy q_ 

i oi 9 

m z jjl z 



(2.40) 
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The structure functions Ffq F are obtained from the F 1 A q YMP via the mass 
factorization relations in Eqs. 



tt __. o-ASYMP j 



( P.10| ) by making the replacements 
iifc -> n^- -, iv ijfe -» ^.- fe -- auu Cf fe SN -> C i>fc on the left and right- 
hand sides respectively. Furthermore we have used the definitions for the 
parton densities in Eqs. ( 2.221) - fl2.24|) . Notice that if the coefficient func- 
tions, indicated by L i k ' , are removed from F F q F one obtains the 
structure functions in the ZM-CSN. The light parton components of the 
structure functions become 



if MSN < LIGHT ( % ,A,Q 2 ,m 2 ) 



fc=i 



III Lb 



1=1 



+f!(n f + l,^®5Cl(n f ,A,^9l) 



+fk+k(nf + 1, V 2 ) ®SC?*(n f , A, 9_, % 



(2-41) 



with 



8C itk (nf,A 



i 9 ) o 



Q' 



C hk (n f + l,^)+C^r(n 



fi' 



„ 91 91 

f ' 9 > 9 

77T yLT 



, rSOFT 



(n /; A. 



0^ tf; 

m 2 ' /z 2 



/fVIRT,ASYMP / Q_ Q_ 



rASYMP/^ V V \ 



(2.42) 



The coefficient functions above satisfy the property that 



Jim 2 5 C, fc (n /; A,— ,— 



J 3>m^ 



/? / , -i Q \ r HARD,ASYMP 



( A Q 2 Q 2 

v m 1 /i 2 



(2.43) 



23 



Using the relations in Eqs. ( |2.9| ), Q2.12J ) one can make a comparison between 



the CSN and the BMSN schemes. In the asymptotic limit the heavy quark 
components satisfy (see Eqs. (^), ( P^| ), ( ggD ) 



lim i^ MSN (n / + l,A,Q 2 ,m 2 ) = lim F^(n f + 1, A, Q\m*) 

Q 2 >m 2 Q 2 ^>m 2 

= lim F™ F (n f + l,A,Q 2 ,m 2 ), (2.44) 

Q 2 ^>rn 2 

and the light parton components satisfy (see Eqs. fl2.30|) , ( |2.41| )) 

lim i^ BMSN ' LIGHT (n / ,A,Q 2 ,m 2 )= lim if SN ' LIGHT (n /; A, Q 2 ,m 2 ) , 

Q 2 ^>m 2 Q 2 ->m 2 

(2.45) 

provided we impose the same boundary conditions on both schemes. From 
the discussion above we infer that the only difference between the CSN and 
the BMSN schemes can be attributed to the m 2 /Q 2 -terms which are present 
in Cj C q N . Such terms do not occur in the C^ q appearing in Ffq ¥ if one chooses 
the BMSN-scheme. In the next Section we will make a study where in Q 2 
these differences are noticeable. 
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3 Comparison between the CSN and the BMSN 
scheme 

In this Section we will make a comparison in the case of charm production 
between the above CSN scheme and the BMSN scheme proposed in ||, [[]]. 
For that purpose we construct a parton density set with four active flavors 
from an existing three flavor set in the literature following Eqs. (2.22)-(2.24). 
The charm quark density of our set will be compared with those in other 
sets with four active flavors presented in |L8| (MRST98, central gluon) and 
[ ITU (CTEQ5HQ). Using our set we will study the differences between the 
charm components of the structure functions Fj C c SN (n/ + 1) in Eq. ( |2.29|) and 
Fj B c MSN (nj + 1) in Eq. ( 2.39| ) in particular in the threshold region. 



Since all coefficient functions are computed in the MS-scheme we choose 
the leading order (LO) and next-to-leading order (NLO) parton density sets 
presented in [[H| which contain three active flavors only (i.e. u,d,s). This 
implies that one has chosen rif = 3 for the anomalous dimensions. In order 
to make this paper self contained we give some details here. In LO where the 
input scale /^o is chosen to be /z 2 , = ^to = 0-26 (GeV/c) 2 the input parton 
densities are 

xu v {3, x, /4 ) = 1.239 x - 48 (1 - x) 2 - 72 (1 - 1.8^ + 9.5x) 

xd v (3,x, n1 ) = 0.614 (1 — x) ' 9 xu v (3,x, hIq) 

xA(3,x,fx 2 LO ) = 0.23 x°- 48 (l-x) 1L3 (l- 12.0^ + 50.9x) 

x(u + d)(3,x,fi 2 LO ) = 1.52 x - 15 (1 - x) 9 ' 1 (1 - 3.6^ + 7.8x) 

xg(3,x,fjl ) = 17.47 x L6 (l-x) 3 - 8 

xs(3,x,fi1 Q ) = xs(x, /i 2 o ) = (3.1) 

In NLO where the input scale equals /i^LO = 0-^0 (GeV/c) 2 we have 

xu„(3,x,/4 LO ) = 0.632 x - 43 (1 - x) 3 - 09 (1 + 18.2x) 

xd v (3, x, /4lo) = 0- 624 (1 - X Y'° xu v (3, x, /4lo) 

xA(3,x,/4 LO ) = 0.20x a43 (l-x) 12 - 4 (l-13.3 v / x + 60.0x) 

x(w + d)(3,x,/4 LO ) = 1.24 x - 20 (1 -x) 8 - 5 (1-2.3^+ 5. 7x) 

xg(3,x,fi 2 NLO ) = 20.80 x L6 (l-x) 4 - 1 

xs(3,x,/4lo) = xs{x,fi^ LO ) = 0. (3.2) 



25 



where A = d — u. Furthermore in |T2j the heavy quark masses are m c = 
1.4 GeV/c 2 , nib = 4.5 GeV/c. In both sets the densities are evolved from 
a very low starting scale, where it is necessary to use the exact numerical 
solution for the running coupling constant a s (/i 2 ). The latter follows from 
the implicit equation 



Q 2 


- 4 " ft In 


4vr A" 
_A)a s (/i 2 ) /?o. 


A MS 


P a 3 {fj, 2 ) ffi 


o = H- 


- \n f , A = 102 


38 
3 '' 



(3.3) 

and will be used in all the following formulae. Furthermore we adopt the 
values A^ 56 = 299.4,246,167.7,67.8 MeV which yield a s (5,M|) = 0.114. 
Notice that the values for A n follow from the matching conditions 

a s (n f ,A nf ,m 2 ) = a s (n f + 1, A nf+1 , m 2 ) . (3.4) 

where for rif = 3 and rif = 4 one has to choose m = m c and m = rrib 
respectively. For the computation of i^ c XACT (gj) and F^ YMP (( ^39l) we 



take rif = 3 for the parton densities and the running coupling constant in 
Eq. (fO[). However for F^ C SN ( pp , i^ B c MSN (|2~39|) and i^ p c DF O) we 



need n/ = 4 for the coupling constant in Eq. ( |5T5| ) and a parton density 
set with four active flavors (i.e. u,d,s,c) when the scale \i becomes larger 
or equal to the heavy flavor mass m. For reasons which will be explained 
below our computations are performed with parton densities represented in 
LO, NLO and NNLO (next-to-next-to leading order). The LO densities 
are convoluted with the order a 2 contributions to the coefficient functions. 
The NLO densities are convoluted with the order a s parts of the coefficient 
functions. The zeroth order contributions to the latter have to be multiplied 
with the NNLO densities. Notice that for nf = 3 we only need LO and 
NLO densities here since the heavy quark coefficient functions in F^ ACT 
and ^asymp s ^ ar ^ j n orc [ er q; s . Our sets with four active flavors are derived 
from the ones with three active flavors by putting n/ = 3 in Eqs. ( p.25|) 



( |2.27| ) starting at a specific scale which we choose as /i 2 = m 2 . Hence it 
follows that in LO or in zeroth order a s one gets 

04,*,m 2 ) = 03,*,m 2 ), 
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/ C L + U e (4, x, m 2 ) = 0, 
f s ' LO (4,x,m 2 ) = f s > LO (3,x,m 2 ), (3.5) 



ig \-1-1 ) Jg 

whereas in NLO or in first order a s one obtains 

/f + L f(4,x,m 2 ) = /JJf(3 ) x,m 8 ), 

/ c N i (4,x,m 2 ) = 0, 

/ s s ' NLO (4,x,m 2 ) = /s s ' NLO (3,x,m 2 ). (3.6) 

Since the two-loop OME's A§i, A^Jg in Eqs. (P^5|)-(P^7D do not vanish at 
/i 2 = m 2 in the MS-scheme we find that in NNLO or in order a 2 the parton 
densities are discontinuous at /x 2 = m 2 while going from three to four flavors 
i.e. 



/fc N + r(4,x,m 2 ) 


± 


/EF°(3,z,m 2 ) 


/™ LO (4,*,m 2 ) 


± 


o, 


/f NLO (4,*,m 2 ) 


± 


ff^°(3,x,m 2 



(3.7) 

Note that if we would drop the terms independent of In fi 2 /m 2 in the two- 
loop operator matrix elements the inequalities in Eq. ( |3.7| ) would become 
equalities. Above ft = m c all four flavor number densities evolve with rif = 4 
(we have not yet included a bottom quark density above \x = m^) . The evolu- 
tion of the parton densities, given by the renormalization group equations in 
Eq. (|2.5| ), is determined by the anomalous dimensions 7},- (LO), 7^ (NLO) 

(2) (2} 

and jlj (NNLO). Unfortunately the three-loop anomalous dimensions 7^ 
are not known yet except for the moments N = 2, 4, 6, 8 (see [p0}1 ). However 



an analysis of the light parton structure function in Eq. ( |2.1|) [pi}] reveals 

(2) 
that the contribution from 7^ is less important numerically than the con- 



tribution due to the two- loop coefficient functions computed in |}13|1 . Hence 
our ignorance about the three-loop anomalous dimension will not apprecia- 
bly alter our results. Therefore our NNLO analysis is only determined by 
the boundary conditions in Eq. ( |3.7| ) which only affects the charm density 
appearing in Fp^ N , i^ B c an d -Ff c DF - The evolution of the parton densities 
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above /i 2 = m 2 , presented in |j22| , was performed using a computer program 
which implements the direct x-space method (similar to that of fl23|). The 
code is written in C++ and has the options to evolve densities in LO and 
NLO whereas the NNLO option presently only uses the NLO anomalous di- 
mensions. We have checked that the evolution of the parton densities is in 



agreement with the results in [24| . 



Before substituting the parton densities into the structure functions we 
encounter a problem caused by the inequalities in Eq. (|3.7|). This happens in 
the threshold region which, according to the ACOT boundary conditions in 
Eq. ( |2.35| ), is defined by Q 2 < m 2 . In this region one has to choose ft 2 = m 2 
so that F^(n f = 4) and if c MSN (n/ = 4) are equal to F^ ACT (n f = 3). 
Notice that the latter has to be understood in the sense mentioned below 
Eq. ( |2.38|) where L,^ — > Lf AKD . However since a s {m 2 ) is rather large 



we have to truncate the perturbation series for the structure functions to 
the desired order otherwise the threshold behavior of all the Fp^ N (nf = 4) 
and _F B c MSN (n/ = 4) will be destroyed. Let us explain this for the BMSN 
scheme in Eq. ( |2.39| ). The arguments for the CSN scheme proceed in an 
analogous way. The conditions that the i^ BMSN (ri/ = 4) = i^ E c XACT (ra / = 3) 
for /i 2 = m 2 implies that the Fff YMP (rif = 3) in Eq. ( j2.39| ) are canceled by 
the FfP F (rif = 4). Further we have to bear in mind that the i^ E c XACT (nj = 3) 
and Fff YMP (rif = 3) are determined by the parton densities in the three 
flavor number scheme whereas the Ff^ ¥ (rif = 4) are constructed out of the 
four flavor number scheme parton densities. The latter have the form 



/ fc (4,x,m 2 ) NNLO = / fc (3,*,m 2 ) NNLO [l + 0(c* 2 



/ 5 s (4,x,m 2 ) NNLO = /g s (3,*,m 2 r LO [l+0(a 2 )] 

c(4,*,m 2 ) NNLO = / s s (3,x,m 2 ) NNL0 0(a 2 )+/ g S (3,x ; m 2 ) NNL0 0(a 2 ). 

(3.8) 



If these densities are substituted in F i p c DF (n/ = 4) in Eq. fl2.40|) one obtains 



additional terms of order a^ and a^ in the perturbation series which are not 
canceled by Fff YMP (nf = 3). Notice that the latter are only computed up 
to order a 2 . This effect is caused by multiplying the four flavor number 
densities with the coefficient functions C^k corrected beyond zeroth order in 
a s . To avoid these higher order terms we propose the following formulae for 
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the structure functions in the CSN scheme 



Ffi»{n f + 1,A,C? 



m 



Q 



f NNLO/ ,-, , 2^CSN,NS,(0)/Q 2 x 

JQ+Q [n f + l,fi )<-i,Q (—5) 



m z 



+a s (n f + l,^)\f^(n f + l,i, 2 )®el 



7 CSN,NS,(1) 



^77l 2 ' /i 2 



+^ na (»/+i,^«c£r*<"(^,^ 



+a>, + 1,/, J ) 0( n/ + 1,a 2 ) ® c, 



CSN,NS,(2)/ . -, Q Q \ 



<' PSl(2) (S|)) + |»/-M^ 2 



C 



CSN 



,PS,(2) / Q Q 

^m 2 ' /x 2 



+^( n/ +i,^«cg»*»(^, 5: 



+aj(n, + 1, /i 2 ) £ 4 /&(», + 1, ^) ® L™' NS ^ (A, £L) . (3.9) 



fe=i 



m" 



Notice that from the perturbative point of view the heavy quark density 
fo+Q starts in order a s (/i 2 ) so that after multiplication with C iQ ' and 

C iQ ' the product becomes of order ai?(/i 2 ). Hence these coefficient func- 

tions did not appear in the mass factorization relations (|2.17 )-( 2.20 ) since 
the latter are carried out up to order a 2 (/i 2 ). Since the heavy quark density 
in CSN has to be treated on the same footing as the light flavor densities, in 
particular after resummation of the terms in In* (/z 2 / m 2 ) , all densities are con- 
sidered to start in zeroth order in perturbation theory. This explains the form 
of of the above expression. Furthermore the coefficient functions C i q ' 

and Ci n ' which originate from parton processes with an heavy quark 

in the initial state have not been calculated yet. Therefore we have to ap- 
proximate them by the replacements 
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respectively. The remaining CSN scheme coefficient functions can be com- 
puted via the relations in Eqs. (|2.17|) - (|2.20| ), which are defined in terms of the 
known iJ's and A's. The light partonic parts of the CNS structure functions 
up to second order are given by 
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If we choose the maximum A = s (defined in the figure caption for Fig. |5|), 
we get L iq ' =0. On the other hand if the minimum value is adopted 

i.e. A = Am 2 one gets L iq ' '*■ ' — L iq and F®q ACT becomes equal 
to the conventional expression given in Eq. (|2.6|) . The structure function 
Ffq YMP is obtained from the expression above by replacing the exact co- 
efficient functions by their asymptotic analogues. Furthermore to calculate 
f bmsn in Eq (p|) we need 
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Finally up to order a 2 , Eq. ( f2.41| ) becomes 
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Notice that we have the relations 
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The first and second terms in the expression above cancel the third last and 



SOFT,NS,(2) 



final term in Eq. ( |3.15| ). The result is then equal to C iq q ' in Eq 



( |3.11|) in the limit Q 2 3> m 2 . The form of the above structure functions 
also suppresses higher order terms beyond a 2 arising from the three flavor 
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number parton densities since the latter also contain terms proportional to 
a s and higher. This becomes apparent if one takes the Nth moments of the 
densities. For instance we observe the following behavior up to NNLO in the 
non-singlet case 
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The choice of the multiplication rules above avoids the appearance of scheme 
dependent terms beyond the order in which we want to compute the structure 
functions. The above prescription guarantees that for Q 2 < m 2 we satisfy 
the condition F^ c = Fp c ; (tj/) in both schemes. 

In the subsequent part of this section we will only discuss the case where 
the heavy quark is the charm quark, i.e. Q — c. Hence in all expressions 
above we have to choose n/ = 3. Further we have to make a choice for 
the cut off A appearing in the coefficient functions L k ' ,( ' (k = 2, L). 

The 5 = (A — 4m 2 )/(s — 4m 2 ) dependence of xL 2>q ' is shown in Fig. 

8, where one notes that it peaks at large x, i.e., near threshold. (The plot 
for xL Lq ' has a similar shape). After convoluting this function with 

the parton densities its contribution to the structure function i^ IGHT only 
amounts to a few percent at Q 2 = 10 3 (GeV/c) 2 . At decreasing Q 2 the 
contribution becomes even smaller. The same holds for L^ RD contributing 
to i*2 c . Hence the dependence of the structure functions on the value of A 
will be very small. Therefore in the subsequent analysis we choose A = 10 
GeV 2 . Other choices hardly affect the plots so that our conclusions will be 
unaltered. 

Next we present the x-dependence of the NNLO charm density (see Eq. 
( |3.7|) ) for various values of /z 2 in Figs. 9a,b. The latter plot emphasizes the 
region 0.01 < x < 1. At \x = m it becomes negative for x < 0.007 which is 
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due to the boundary condition in Eq. ( f2.23j ) and the momentum sum rule. 



When /j, > m the density becomes positive over the whole x range. In Figs. 
9c and 9d we have shown the charm densities in NLO which are obtained 
from [[18] and |19] respectively, with an offset scale so that they can easily be 



compared with Fig. 9a. The former is constructed in the TR scheme whereas 
the latter follows the prescription of ACOT. Both are positive over the whole 
x range. Our LO and NLO parameterizations, which are not shown in the 
figures, are also positive for all values of x. This property can be traced 
back to the boundary conditions which yield in LO and NLO c(x,m 2 ) = 0. 
Note that a direct comparison between the charm densities from different 
groups is not meaningful because each group fits different data to determine 
their respective input three flavor number gluon densities. However it seems 
that our charm density, at small /i 2 , does not rise as steeply as that of the 
CTEQ5HQ [|Tj| at small x. It is more similar to the MRST98 (central gluon) 
[ |T8|j density. In Fig. 10 we make a comparison between our charm density 
which evolves according to the renormalization group equation and the one 
computed in fixed order perturbation theory (FOPT) via Eq. ( p.23j ) . To that 
order we have plotted 

c EVOLVED/ x 2\ 

R(x, /i ) = ^pt/ 2 x • (3-19) 

The density c FOPT is computed up to order a 2 since the OME's in Eq. ( J2.26Q 
are only known up to that order. In Fig. 10a we have shown the ratio in LO. 
The latter implies that we have only kept the leading logarithms in In /i 2 /m 2 
in the OME's which are resummed in all orders in c EVOLVED . The deviation of 
R from unity shows the effect of the resummation. The same ratio is shown in 
NLO in Fig. 10b where we also included the subleading terms in the OME's. 
Finally if we take into account the non-logarithmic terms in the two-loop 
OME's A^f and Aq S / 2) ( ^26|) one obtains the NNLO ratio (see Fig. 10c) 



The figures reveal that in LO and NLO the effect of the resummation is very 
small except near x — 1. This picture changes if we go to NNLO where the 
deviation of R from one becomes appreciable when x tends to zero. Here 
R can even become negative which happens for /i 2 p» 3 (GeV/c) 2 . This 
effect is wholly due to the boundary condition c(x, m 2 ) ^ which occurs 
beyond NLO. Furthermore the figures reveal that R > 1 at large x whereas 
R < 1 at small x. Notice that in Fig. 10c c FOPT (x,/i 2 ) = for x = 0.007 
at /j, 2 = 2 (GeV/c) 2 so that R = oo which explains the bump in the figure. 
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Figure 10c is important because it shows that c E VED (x, fj, 2 ) < c PT (x, /i 2 ) 
at small x. The consequence is that F B c MSN (x, Q 2 ) and Fj C c SN (x, Q 2 ) will 
become smaller than F E c XACT (x, Q 2 ) when Q 2 becomes slightly larger than 
m 2 due to the choice made for the factorization scale in Eq. ( |2.35| ). This can 
even lead to a negative structure function as will happen for F^ N which we 
will see later on. 

Now we present results for the various structure functions. In Fig. 11 we 
show the charm quark structure functions in NNLO given by F^^i^f = 4), 
f bm SN( n/ = 4 ) 5 F 2 p c DF (n/ = 4) and F^ ACT (n f = 3) plotted in the region 
1 < Q 2 < 10 3 (GeV/c) 2 for x = 0.05. The figure reveals that there is 
hardly any difference between the BMSN and CSN prescriptions. The curves 
in both prescriptions lie between the ones representing i r 2 > c )F ( n / = 4) and 
F 2 EXACT (?2j = 3) except for low Q 2 . In this region the latter is a little 
bit larger than the other ones which is expected from the discussion of the 
charm density given above. Notice that in the low Q 2 region F 2 P EF (^/ = 4) 
becomes negative which means that charm quark electroproduction cannot 
be described by this quantity anymore. In Fig. 12 we present the same 
plots for x = 0.005. Again one cannot distinguish between F 2 B ( ^ ISN (n/ = 4) 
and i 7 2 C f N ( ?7 '/ = 4) but now both are smaller than F^ ( n f = 3) over the 
whole Q 2 range. The latter is even larger than F p ^ (nf = 4) in particular 
for Q 2 > 50 (GeV/c) 2 . Further we want to emphasize that due to our 
careful treatment of the threshold region there is an excellent cancellation 
(to three significant places) between i 7 2 > <? F ( ri / = 4) and and F A ^ YMP (nf = 3) 
at very small Q 2 so that both F^ N (n f = 4) and F 2 B c MSN (n/= 4) tend to 
F EXACT (nj = 3). Also at large Q 2 we have an excellent cancellation between 
f Asymp^ / = 3 ) and F EXACT( n/ = 3) so that both F™ N (n f = 4) and 

F 2 B c MSN (n/ = 4) tend to F 2 p c DF (n/ = 4) (see Eq. (ggj) ). 



In Fig. 13 we show similar plots as in Fig. 11 but now for the charm 
quark longitudinal structure functions. Here we observe a difference between 
the plots for F^(n f = 4) and F^ SN (n f = 4) in the region m 2 c < Q 2 < 
40 (GeV/c) 2 . In particular the latter tends to F^ ACT (n f = 3) while the 
former is larger. Furthermore F PBF (n/ = 3) is considerably larger than 
the other three structure functions, which differs from the behavior seen in 
Fig. 11. This can be mainly attributed to the gluon density which plays 
a more prominant role in Fi jC than in F 2iC . For x = 0.005 (see Fig. 14) 
the difference between the BMSN and the CSN descriptions becomes even 
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more conspicuous. In this case F^ N (nf = 4) becomes negative in the region 
m 2 < Q 2 < 7 (GeV/c) 2 which is unphysical. This effect can be attributed to 
the zeroth order longitudinal coefficient function in Eq. ( |3.9p , which behaves 
like C L q ' — Am 2 /Q 2 (see [0), and the non-vanishing charm density at 

/i 2 = m 2 . In the case of BMSN the longitudinal coefficient function is equal 
to zero in lowest order so that F^^ (rif = 4) does not become negative. 

In Figs. 15 and 16 we make a comparison between the NLO and the 
NNLO structure functions F 2 c c SN (n/ = 4) and F 2 B c MSN ( n / = 4 )- Both P re " 
scriptions i.e. CSN and BMSN lead to the same result in NNLO. However 
while going from NLO to NNLO the the structure function F 2( ? (n/ = 4) 
decreases whereas F^^ i n f = 4) increases a little bit. The differences in 
the case of x — 0.005 are even smaller than those observed for x = 0.05. 
The same comparison between NLO and NNLO is made for the longitudi- 
nal structure functions in Figs. 17 and 18. Here the differences between 
NLO and NNLO are much larger than in the case of F 2jC in Figs. 15,16. In 
NLO Fl^ SN (n f = 4) is smaller than the one plotted for NNLO. However 
for F^ N (nf = 4) we see a decrease in the small Q 2 -region while going from 
NLO to NNLO whereas for large Q 2 we observe the opposite. In particular 
the valley in the region m 2 < Q 2 < 7 (GeV/c) 2 observed for F^ K (nf = 4) at 
x = 0.005 in NNLO turns into a bump. This is due the boundary condition 
on the charm density which in NLO vanishes at /i = m c whereas in NNLO 
it is negative at small x- values (see Fig. 10c). From the observations above 
one can conclude that the CSN prescription bedevils the threshold (low Q 2 ) 
behavior for Fl <c due to the non-vanishing zeroth order longitudinal coeffi- 
cient function C L q ' . This problem is avoided by TR in [10] by imposing 



a condition on the structure functions as indicated in Eq. ( |2.36| ). Hence 
our results for i^ B c agree reasonably well for i = 2 and i = L with those 
presented in NLO by TR in ||10| . This is mainly due to the fact that there is 



only a small difference between the NLO and NNLO approximations in the 
BMSN scheme. It also reveals that the condition in Eq. ( 2.36|) for i = L can 



be mimicked by a vanishing zeroth order longitudinal coefficient function. 
Note that results for the x-values presented above are representative for the 
whole range 5 x 10~ 5 < x < 0.5. 

To summarize the main points of this paper we have discussed two vari- 
able flavor number schemes for charm quark electroproduction in NNLO. 
They are distinguished by the way mass factorization is implemented. In the 
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CSN scheme this is done with respect to the full heavy and light quark struc- 
ture functions at finite Q 2 . In the BMSN scheme the mass factorization is 
only applied to the coefficient functions in the large Q 2 limit. Both schemes 
require three flavor and four flavor number parton densities which satisfy 
NNLO matching conditions at a scale /i 2 = m 2 . We have constructed these 
densities using our own evolution code. The schemes also require matching 
conditions on the coefficient functions which are implemented in this paper. 
We have also made a careful analysis of the removal of dangerous terms in 
ln(Q 2 /m 2 ) from the Compton contributions so that both Fp^ N (rif + 1) and 
F i B c MSN (rij + 1) are collinearly safe. We have done this in a way which is 
simplest from the theoretical point of view, by implementing a cut A on 
the mass of the c — c pair which has to be determined by experiment. This 
cut is not required in the fixed order perturbation theory approach given by 
i^ E c XACT in B for moderate Q 2 -values. 

Finally we made a careful analysis of the threshold behaviors of Fp^ N (nf+ 
1) and -Fj B c MSN (ri/ + 1). In order to achieve the required cancellations near 
threshold so that they both become equal to F^ ACT (rif) one must be very 
careful to combine terms with the same order in the expansion in the running 
coupling constant a s . Therefore technically we require six sets of parton 
densities, namely the LO, NLO and NNLO three flavor number sets and the 
LO, NLO and NNLO four flavor number sets. However not all the necessary 
theoretical inputs are available to us to finish this task. The approximations 
we made in this paper were sufficient to provide very clear answers. We 
successfully implemented the required cancellations near threshold and the 
corresponding limits at large scales came out naturally. Inconsistent sets of 
parton densities automatically spoil these cancelations. We did not have to 
use matching conditions on derivatives of structure functions as proposed 
in |TIJ, which seem very artificial. The numerical results do however end 



up quite similar. We have also shown that the CSN scheme defined above 
leads to an unnatural behavior of the longitudinal structure function in the 
threshold region which is due to a non-vanishing zeroth order coefficient 
function. Since there are no other differences between the CSN and BMSN 
schemes we recommend the latter because it is less complicated than the 
former. In particular it does not need additional coefficient functions other 
than the existing heavy quark and light parton coefficient functions available 
in the literature. 
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Appendix A 



In this appendix we present the exact expressions for the heavy quark coef- 
ficient functions L\ corresponding to the Compton process in Fig. |5| when 
there is a cut A on the invariant mass Sqq of the heavy quark pair. As 
explained in [I5j the calculation is straightforward because one can first inte- 
grate over the heavy quark momenta in the final state without affecting the 
momentum of the remaining light quark. The phase space integrals are the 
same as the ones obtained for the process 7(g) + q(k 1 ) — > g* + q (g* — > Q + Q) 
where the gluon becomes virtual. In the expressions for the complete inte- 
gration over the virtual mass Sqq of the gluon (see Fig. |) one integrates 
over the range 4m 2 < Sqq < s with s = (q + kif. The resulting expressions, 



called Li with r = NS, S and i = 2, L are presented in appendix A of 
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where the partonic scaling variable is equal to z — Q 2 /(2q- ki) = Q 2 /(s + Q 2 ), 
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The variable A, which allows us to distinguish between soft and hard (ob- 
servable) heavy quark anti-quark pairs, is in the range 4m 2 < A < s. The 
variable z is in the range < z < £/(£ + 4). 

Note that when A = s one obtains L Lq ' — > L\^ q which are re- 



ported in p5| . When the integration range is given by A < Sqq < s we get 
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Notice that L, „ 'is finite in the limit m — > so that it does not contain 
collinear divergences. The latter can be wholly attributed to L iq ' ' ( ' as 
is revealed if one takes the limit Q 2 —>■ oo. In this case the expressions (A.l) 
and (A. 2) reduce to 
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respectively. In the limit m — > the results above show the same logarithmic 
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terms in W(Q /m ) (i = 1, 2) as the asymptotic expressions for Lj'~ given in 
Eqs. (D.7) and (D.8) of [15]. Hence the differences between the results there 
and the asymptotic expressions above are free of any terms in ln(Q 2 /ra 2 ). 
The asymptotic expressions of L { ' for Q 2 — > 00 are given by 
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respectively. As has been mentioned below Eq. (|3|) the expressions above 
are finite in the limit m — ► (b — > 1) so that they do not contain collinear 
divergences. 
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Figure Captions 
Fig. 1. The lowest-order photon-gluon fusion process 7* + g — > Q + Q con- 



tributing to the coefficient functions H, 



s,(i) 



Fig. 2. Some virtual gluon corrections to the process 7* + g — »■ Q + Q 
contributing to the coefficient functions if^ . 

Fig. 3. The bremsstrahlung process 7* + g — *■ Q + Q + g contributing to the 



coefficient functions H , „ \ 

Fig. 4. The Bethe-Heitler process 7* + q(q) ^ Q + Q + q(q) contributing 
to the coefficient functions H iq . The light quarks q and the heavy 
quarks Q are indicated by dashed and solid lines respectively. 

Fig. 5. The Compton process 7* + q(q) — ► Q + Q + q(q) contributing to the 
coefficient functions L iq . The light quarks q and the heavy quarks 
Q are indicated by dashed and solid lines respectively (s = (p + q) 2 , 
Sqq = (Pi +P2) 2 see text). 

Fig. 6. The two-loop vertex correction to the process 7* + q — > g contain- 
ing a heavy quark (Q) loop. It contributes to C iq ' (Q 2 /m 2 ) = 



F^{Q 2 /m 2 ) Cg> 






Fig. 7. Order a s corrections to the process 7* + Q — > Q an< i the reaction 
7* + Q — > <5 + 5 1 contributing to the coefficient functions H iQ ,( ' . 

Fig. 8. The <5 = (A— 4m 2 )/(s— 4m 2 ) dependence of xL 2) q ' (x,Q 2 /m 2 ,A) 
at Q 2 /m 2 = 50 (Eq. (A. 2)) plotted as a function of x for <5 = 1, 0.1, 
0.01 and 0.001 respectively. 

Fig. 9. (a) The charm density xc NN (4, x, /x 2 ) shown in the range 10 -5 < 
x < 1 for fi 2 = 1.96, 2, 3, 4, 5, 10 and 100 in units of (GeV/c) 2 . (b) 
similar plot as in (a) but now for 0.01 < x < 1. For a comparison we 
have also shown the NLO results obtained by MRST98 and CTEQ5HQ 
for the range 10~ 5 < x < 1 in (c) and (d) respectively. 

Fig. 10. Ratios R(x, /1 2 ) = xc EVOLVED (4, x, /i 2 )/xc FOPT (4, x, /1 2 ) for the scales 
/1 2 =2, 3, 4, 5, 10, 100 in units of (GeV/c) 2 . (a) LO, (b) NLO, (c) 
NNLO. 
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Fig. 11. The charm quark structure functions F^ T (rif = 3) (solid line) 
F 2 c c SN (n/ = 4), (dot-dashed line) F 2 B c MSN (n/' = 4), (dashed line) and 
F 2 p ( P F (nj = 4), (dotted line) in NNLO for x = 0.05 plotted as functions 

of'g 2 . 

Fig. 12. Same as in Fig. 11 but now for x = 0.005. 

Fig. 13. The charm quark structure functions F^ T (n/ = 3) (solid line) 
Fg N (n/ = 4), (dot-dashed line) F^ SN (n f '= 4), (dashed line) and 
Ff^irif = 4), (dotted line) in NNLO for x = 0.05 plotted as functions 

ofg 2 . 

Fig. 14. Same as in Fig. 13 but now for x = 0.005. 

Fig. 15. The charm quark structure functions F™ (rif = 4) in NLO (solid 
line), NNLO (dotted line) for x = 0.05 and F 2 c c SN (n/ = 4) in NLO 
(dashed line), NNLO (dot-dashed line) for x = 0.05 plotted as functions 

ofg 2 . 

Fig. 16. Same as in Fig. 15 but now for x = 0.005. 

Fig. 17. The charm quark structure functions F^^ ISN (n/ = 4) in NLO (solid 
line), NNLO (dotted line) for x = 0.005 and F^ N (n f = 4) in NLO 
(dashed line), NNLO (dot-dashed line) for x = 0.05 plotted as functions 
of Q 2 . 

Fig. 18. Same as in Fig. 17 but now for x = 0.005. 
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